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AN APPROACH- GUIDANCE METHOD USING A SINGLE 
ONBOARD OPTICAL MEASUREMENT 


By Harold A. Hamer and Katherine G. Johnson 
Langley Research Center 

SUMMARY 

An empirical method is developed for onboard guidance within the sphere of influ- 
ence of a celestial body. The procedure requires only limited onboard calculations and 
leads to approach-guidance predictions sufficiently accurate for emergency or backup 
operations. The method is applied to lunar-approach trajectories and is studied in detail 
for certain lunar missions. 

The procedure relies heavily on use of precalculated data and is unique in that only 
a single angular measurement from the star to the moon is required, provided that it is 
made at or near the lunar sphere of influence. If the approach guidance is delayed to a 
time closer to the moon, an additional measurement of the subtended angle of the moon is 
required. The method is designed specifically to control the magnitude of the perilune 
radius, but the perilune position and velocity values automatically remain close to the 
nominal values. An error analysis with an assumed one-sigma error of 10 seconds of arc 
in the optical angular measurements and a one-sigma velocity-cutoff error of 0.2 m/sec 
has shown that the perilune radius can be controlled to a one-sigma accuracy of from 7 to 
13 km, depending on the time the approach guidance is performed. The effects of maneu- 
vering errors, star location, and empirical approximation errors on the approach guidance 
are discussed. 

For manned flight, the required angular measurements can be readily made from 
onboard the spacecraft by a sextant -type instrument. For unmanned flight, the measure- 
ments can be made automatically by pointing the spacecraft (or tracker instrument) in a 
predetermined direction to a star and then sighting to the planet with a scanner-type 
instrument. 


INTRODUCTION 

In space missions the navigation and guidance is normally accomplished by auto- 
matic procedures which employ earth-based radar measurements. The inclusion of 



procedures which utilize onboard measurements is a desirable feature, both for manned 
and unmanned flights. For example, in interplanetary flight, the planet ephemeris error 
can lead to unacceptable errors in the trajectory position relative to the planet when only 
earth-based measurements are used. This position error is most significant during the 
approach phase, that is, when the spacecraft is within the sphere of influence of the planet. 
Some type of onboard measurements relative to the planet may be required to correct this 
position error. 

Over the years a number of studies have been made to develop onboard guidance pro- 
cedures for controlling the approach to a celestial body. (For example, see refs. 1 to 6.) 

In general, these methods require several types of measurements and the measurements 
must be repeated a number of times. The guidance correction is ordinarily based on a 
statistical filtering technique which requires extensive calculations and computing equip- 
ment, and the methods may also require more than one guidance maneuver. 

The method presented in this paper is unique in that only a single onboard position 
fix is required to determine the guidance correction. This fix is made at a chosen time 
and, when used in conjunction with some simple empirical approximations derived from 
two-body theory, is sufficient to control the periapsis magnitude with a reasonable degree 
of accuracy by using only one guidance maneuver. In this paper, the method is studied in 
regard to the approach phase of earth-moon trajectories, but it can also be applied to 
reentry control for moon-earth trajectories. (See ref. 7.) The method warrants investi- 
gation for application to interplanetary-trajectory control. 

For guidance at or near the lunar sphere of influence, the only measurement 
required is the included angle between a star and the target body. If the approach guidance 
is made closer to the moon, a subtended-angle measurement is also required for deter- 
mining the range. An analysis is included which shows that for acceptable accuracy, the 
star must be in a specified direction with respect to the nominal trajectory. The direction 
of this star and that of the guidance-velocity vector can be predetermined as can most of 
the calculations for deriving the magnitude of the guidance correction. These character- 
istics make the method particularly useful for manual operation. 

The accuracy characteristics of the method are examined by means of a Monte Carlo 
error analysis. The analysis includes the effects of measurement error, thrust-cutoff 
error, and approximation error which is caused by assumptions in the empirical proce- 
dure. The results were obtained by use of the Jet Propulsion Laboratory n-body trajec- 
tory program. (See ref. 8.) 
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SYMBOLS 


D position deviation in direction of specified star 

i,m,n direction cosine of line of sight to star with respect to X-, Y-, and Z-axis, 
respectively 

R lunar radius 

r range to moon center 

r e range to earth center 

rp perilune radius 

Ar mc incremental range to earth center at time of midcourse-position fix, 
r e,a “ r e,n 

Arp incremental perilune radius, *p,a “ r p,n 

T time from earth injection 

Tp time to nominal perilune time 

Tpf time of midcourse-position fix 

At time increment (appendix A) 

V vehicle velocity 

Vp perilune velocity 

AV guidance -velocity correction 

x,y,z position coordinates in Cartesian axis system in which X-axis is in direction 

of Aries, XY-plane is parallel to earth equatorial plane, and Z-axis is in 
direction of north celestial pole 

x,y,z velocity coordinates in Cartesian axis system 
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position coordinates in rotating Cartesian axis system in which X r -axis lies 
along earth-moon line, X r Y r -plane is in earth-moon plane, and Z r -axis is 
in northerly direction 

semisubtended angle of moon 

angle between nominal range vector and major axis of position-error ellipsoid 
flight -path angle 

angle formed at vehicle between line to star and its projection in the instan- 
taneous earth-moon-vehicle plane 

eccentricity of orbit (appendix C) 

angle between nominal velocity vector and major axis of velocity-error 
ellipsoid 

included angle between star and celestial -body center 
angular measurements used for determining a and 9 

angle formed at vehicle between line to moon center and projection of line to 
star in the instantaneous earth-moon-vehicle plane 

angle between approach-guidance-velocity vector and vehicle-velocity vectoj 

product of universal gravitational constant and mass of moon 

standard deviation or root-moon-square value 

true anomaly 

in-plane angle between second midcourse AVg and V n 
absolute value 
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Subscripts: 


A 

a 

add 

D 

F 

m 

n 

R 

r 

r,mc 

r,P 

S 

s 

u 

V,p 


value immediately following approach-guidance correction 
actual 

additional approach-guidance velocity to account for second midcourse 
manuever 

position deviation 

first midcourse correction 

measured value 

nominal value 

moon radius 

range to moon center 

range to earth center at time of midcourse-position fix 

perilune radius 

second midcourse correction 


position deviation from nominal trajectory, 

— 1/2 
(x a - x n ) 2 + (y a - yn) 2 + (z a - Zn) 2 


velocity deviation from nominal trajectory, 

• ■ 2 . o . . 2 

(x a - x n r + (y a - y n ) + (z a - z n ) 


perilune velocity 
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a 


semisubtended angle of moon 


9 star-to-body angle 

A bar over a symbol indicates a vector. 

BASIC METHOD 
Synopsis 

The approach-guidance procedure presented herein is designed to correct perilune- 
magnitude error. The direction of the approach -guidance velocity correction is preset; 
the manuever can be applied at any preselected time within the lunar sphere of influence. 
The direction in which the correction is usually applied is perpendicular to the nominal 
velocity vector and in the orbital plane. This direction is essentially optimum for most 
distances from the moon. 

Midcourse guidance . - The approach-guidance procedure is an outgrowth of the 
onboard midcourse-guidance procedure developed in reference 9 and is designed to 
correct errors resulting from the use of this type of midcourse guidance. In practice, 
the midcourse guidance would be required to correct the actual trajectory in order to 
remove perturbations due to errors attributed to injection and other sources. To reduce 
the amount of analysis required for this paper only one such perturbed trajectory was 
employed for most of the approach-guidance analysis. Some data, however, are given to 
show that the approach-guidance method is capable of handling a family of trajectories 
originating from a wide range of injection errors. The injection errors considered were 
essentially spherically distributed and had one-sigma values of roughly 3 km in position 
and 3 m/sec in velocity. These perturbations are relatively large with respect to 
present-day values. (For the smaller injection errors associated with normal operation, 
the onboard approach guidance could be used in lieu of any other type of midcourse 
guidance.) 

A fixed-time-of-arrival law was used for the onboard midcourse guidance. (See 
ref. 9.) The aim point was selected at the lunar sphere of influence. The first midcourse 
maneuver was performed at 10 hours from injection for 70-hour translunar trajectories 
and at 15 hours for 90-hour translunar trajectories, the position fix being taken 1/2 hour 
before the maneuver. These trajectories span the range of reasonable trip times. 

Because the first midcourse maneuver corrects only for the position error at the aim 
point, a second midcourse maneuver is normally required at the aim point to correct the 
spacecraft-velocity vector back to the nominal vector. 



Approach guidance .- The approach-guidance method is empirical in nature and most 
of the calculations are preflight calculations. In these calculations the guidance -velocity 
requirements are developed from two-body relationships, wherein a closed-form expres- 
sion relating perilune and upstream conditions can be written. These conditions can be 
related to deviations from the nominal trajectory. The deviations can be determined by 
simple onboard optical angular measurements. 

Results obtained from Monte Carlo samples of trajectories perturbed at first mid- 
course are used to show that within the sphere of influence, the position deviation in a 
certain direction predicts the perilune radius and perilune velocity with relatively high 
accuracy. The guidance velocity required to correct the perilune radius is then deter- 
mined empirically as a function of this deviation. It is this precalculated variation which 
the navigator employs for the onboard approach guidance. The only onboard calculation 
required is the simple computation of the deviation by use of the measured and nominal 
values of the onboard measurements. 

Errors considered in guidance procedures . - The random perturbations after first 
midcourse were assumed to be caused by onboard measurement error alone; the effect of 
maneuvering errors would be negligible. At second midcourse the effect of measurement 
error was considered, inasmuch as this maneuver is derived from the first midcourse 
measurements. Although small, the effect of velocity-cutoff error was included in the 
second midcourse correction. Because of the relatively small second midcourse correc- 
tion required, the effect of error in the pointing direction was neglected. In the approach- 
guidance procedure, types of errors considered were: measurement errors, velocity- 
cutoff errors, and approximation errors caused by use of the empirical procedure. The 
effect of approximation error in the midcourse guidance is essentially eliminated by 
the approach guidance. It should be noted that n-body trajectories are used throughout 
this paper in developing the method and in performing the error analyses. 

General Considerations 

The onboard midcourse-guidance procedure of reference 9 reduces the error at the 
aim point significantly; however, the remaining errors must be reduced further by means 
of some type of approach guidance. The approach-guidance procedure presented herein 
proposes to correct such errors resulting from the use of this type of midcourse guidance. 
Some typical errors resulting from the midcourse procedure are shown in figure 1. The 
data are shown for two magnitudes of range -measurement error inasmuch as this error 
has a predominant effect on the aim -point accuracy. The range of aim points covered 
includes all points from the sphere of influence to perilune. 

The accuracy characteristics of the approach-guidance procedure are dependent 
upon the errors associated with the midcourse procedure. As previously stated, the 
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maneuvering errors at midcourse are negligible compared with the measurement errors 
made in estimating the required midcourse correction. The measurements include a 
range determination and three star-to-body angles. Details of the midcourse procedure 
can be found in reference 9. The facts pertinent to the present study are that the predom- 
inant measurement errors and the position-determination errors are generally in the 
direction of the vehicle-velocity vector. Hence, the midcourse-guidance error is essen- 
tially an error in the magnitude of AV rather than in its direction and can be accurately 
controlled in the approach-guidance procedure. As an example of the midcourse-guidance 
errors, one-sigma values along the three axes of the error ellipsoid relating to the 
covariance matrix of midcourse-velocity errors are 1.38, 0.505, and 0.226 m/sec. 

These magnitudes pertain to the data shown in figure 1 for range-measurement error 
°r,mc = 22 km. The values signify that the error at first midcourse is predominantly 
along the major axis of the ellipsoid; the angle between the major axis and the spacecraft- 
velocity vector is about 2° in this case. 

Thrust Assumptions 

In the approach-guidance procedure, the thrust is considered to be impulsive in 
effect; that is, the burning time is negligible relative to the trajectory time scale. The 
impulsive correction is assumed to be applied in a constant direction in the nominal plane 
of motion at initialization of the thrust maneuver. Unless otherwise noted, the velocity- 
correction vector is perpendicular to the nominal velocity vector. Except for the effect 
of engine-cutoff error, the guidance correction is assumed to be perfectly executed. 

These assumptions are all appropriate inasmuch as their effect on the overall results is 
negligible. Finally, a high-thrust device for implementing the approach-guidance maneu- 
ver is assumed. 


Approach- Guidance Procedure 

The approach guidance can be applied at a predetermined time anywhere within the 
lunar sphere of influence. Most of the required calculations are performed before the 
flight by using information on the nominal trajectory. In regard to the midcourse proce- 
dure, the approach guidance must be applied at or beyond the aim point where the second 
midcourse maneuver is made. These two maneuvers can be conveniently combined even 
though in practice some time is required to make the simple approach-guidance measure- 
ments and calculations. A delay of several minutes in the guidance does not appreciably 
affect the overall accuracy of the system. Longer delays, if necessary, can easily be 
taken into account in the procedure. 
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Measurement equations .- The approach-guidance procedure is based on the deter- 
mination at a given time of the deviation of the trajectory in position from the nominal 
trajectory. (See fig. 2.) As will be shown, for a strong correlation with the perilune 
conditions rp and Vp, the deviation must be measured in a predetermined direction to 
a specified star. From the figure, the deviation D is given by 

D = r m cos 0 m - r n cos 9 n (1) 

where r n and 0 n are nominal values. Normally, two measurements are required: 
the range to the moon r m and the angle 0 m included between the star and the center 
of the moon. In this paper, the range is assumed to be measured by the semisubtended 
angle a (that is, r = R/sin a where R is the known radius of the moon). The two 
measurements must be referenced to a common time. Since it may not be possible to 
make simultaneous measurements, a method for updating the measurements is given in 
appendix A. As will be pointed out, 0 n at the sphere of influence should be approxi- 
mately 90°, in which case the effect of error in range is negligible and 

D w r n ^cos 0 m - cos (2) 

Thus no range measurement is required. 

Guidance -velocity determination .- As previously stated, there is a strong correla- 
tion between the deviation D, taken in a certain direction, and the perilune conditions 
Tp and Vp. This correlation is shown in figure 3. In this figure, D is determined at 
a time near the sphere of influence and the measurement star is in the nominal orbital 
plane and perpendicular to the nominal range vector to the moon. The trajectories A 
and B are perturbed differently at injection; trajectory A resulted in a Ar mc of 
-419 km at the time of the first midcourse-position fix, whereas trajectory B resulted 
in a Ar mc of -1287 km. In both cases, the velocity-error vector was about 17.5° from 
the vehicle-velocity vector; for trajectory B the position-error vector was about 79.5° 
from the vehicle-position vector. Each data point in figure 3 represents the result caused 
by random measurement errors in correcting the trajectories at midcourse. Trajectory 
results are shown for two magnitudes of range-measurement error. The smoothness of 
the data for the two trajectories indicates that the approach-guidance accuracy would be 
insensitive to the magnitude of injection error. 

The strong correlation between D and the perilune conditions is a clue that D 
can be used empirically to determine the guidance-velocity correction required to attain 
the desired perilune radius. The equation which relates the approach-guidance-velocity 
magnitude to the conditions at perilune, and hence to the deviation D is as follows: 
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AV = 


V r 2 cos y cos (y + X) - r ppi 2 cos X 


r p ,n 2 ~ r 2 cos 2 (y + X) 

rp,nf(r 2 - r p ,n 2 ) V 2 sin 2 X + r 2 cos 2 (y + X) - r pjI1 2 


2 ii 2^\W 2 
r p,n" r JJ 


r p?n 2 - r 2 cos 2 (y + X) 


(3) 


where 


cos y = 


r p V P 

rV 


(4) 


The equation was derived from the principle of conservation of angular momentum, as 
shown in appendix B. This equation is not used by the navigator onboard the spacecraft, 
but it is used in a preflight analysis to determine the variation of AV with D. 

In the equation, X is the angle between the AV vector and the vehicle-velocity 
vector, which is a predetermined value; y is the flight-path angle. The quantities r 
and V are the values of range and velocity for the perturbed trajectory at the time of 
the approach correction. The quantity r pn is the nominal or desired perilune radius, 
whereas r p and V P are the values of perilune radius and perilune velocity of the per- 
turbed trajectory. The values of r P)I1 and X are constant. For a given time, changes 
in r and V do not significantly affect AV; the variables r p and V p are the main 
contributors of AV. These variables, in turn, are dependent upon y as is seen in 
equation (4). 

Star direction. - For the approach-guidance method the star must be in a given 
direction with respect to the nominal trajectory, and it is not necessary that the variations 
such as those shown in figure 3 be linear. The approach- guidance accuracy is closely 
related to the amount of scatter in the perilune -radius predictions; hence, it is imperative 
that the scatter be kept to a minimum. Examples of scatter for various directions of the 
deviation D at T D = 14.617 hours are shown in figure 4. It is obvious that minimum 
scatter is obtained for a deviation (star measurement) in or near the nominal orbital plane 
and approximately perpendicular to the nominal range vector. The instantaneous nominal 
earth-moon-vehicle plane is used as the reference plane since it is within 0.667° of the 
selenocentric orbital plane of the spacecraft. 

Figure 5 presents data for various times within the sphere of influence. (Note the 
staggered vertical scale. The scale should be read so that in all cases the curve at 
D = 0 would pass through the nominal value of r p , which is roughly 3404 km.) In 
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figure 5(a) examples are shown where the nominal 9 is held constant; that is, the 
direction of the star is always perpendicular to the nominal range vector at the corre- 
sponding time Tp. The direction of the star changes from case to case, and the scatter 
becomes unacceptable as the approach-guidance measurement is delayed to times near 
the moon. The characteristics of the scatter in figure 5(b) show that the optimum direc- 
tion of the star does not change with time. The indicated change in the true anomaly 
depicts the angle through which the range vector rotates and hence the change in the 
nominal value of 9 with time. For example, at Tp = 14.617 hours, 0 n would be 90°; 
whereas at T p = 4.617 hours, 0 n would be about 80° (or 100°). It is apparent from 
figure 5(b) that regardless of the time that the measurement is made, the star must be 
in a direction perpendicular to the nominal range vector at the sphere of influence. This 
direction is illustrated in figure 2. 

The optimum direction of the star is apparently in the orbital plane and perpendic- 
ular to the nominal range vector at the sphere of influence. The characteristics of the 
scatter in figure 6 show that the star can be as much as 2° away from the optimum direc- 
tion and still give adequate accuracy. Although it is not shown, an in-plane displacement 
from the perpendicular direction much greater than 2° would not be acceptable. Note the 
amount of scatter in figure 4(c) for which the displacement is about 9.5°. Out-of -plane 
displacement up to at least 30° does not appreciably affect the accuracy insofar as scatter 
is concerned; however, it does affect the measurement sensitivity as shown in figure 7. 

From the sketch in the upper part of figure 7, it can be determined from the right 
spherical trigonometric relationship 

cos 9 = cos 0 cos 6 


that 


d 9 _ f cos^5 - cos2g \ ^ 
d® \ 1 - cos^0 / 


where — ^ could be called the measurement- sensitivity factor. For a star direction 
at) 

perpendicular to the orbital plane, this factor is zero; this statement means that a mea- 
surement in this direction gives no indication of the in-plane trajectory deviation, which 
is essential for controlling the magnitude of r p . (The out-of-plane deviation would be 
important for methods which control the location of r p .) Inasmuch as figure 5 shows 
that the nominal value of 9 changes with time T p , the data in figure 7 are shown for a 
change of ±10°. This amount of change shows little or no effect on the sensitivity factor, 
It is seen in figure 7 that an out-of-plane displacement of 30° has only about 10-percent 
effect on the sensitivity. 
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Optimum direc tion and time of AV.- For guidance made within the lunar sphere of 
influence, the optimum direction of the AV vector is in the orbital plane and essentially 
perpendicular to the nominal velocity vector. This condition is illustrated in figure 8 
which is an example of the correction of one perturbed trajectory which had a position 
error of 271 km at the sphere of influence {Tp = 14.617 hr). Note that close to the moon 
(Tp = 2.617 hr), the optimum direction for AV is about 80°; however, the change in AV 
from that at A = 90° is negligible. Also, as expected, the figure shows that for minimum 
fuel requirements, the optimum time for the guidance maneuver is at the sphere of influ- 
ence. The increase in fuel requirements as the maneuver is delayed to times closer to 
the moon is also illustrated in figure 9. The variation shown applies to the energy level 
of the translunar trajectory shown in figure 10 and is only approximate because it was 
determined for one particular perturbed trajectory. However, for any other perturbed 
trajectory, the values would not differ appreciably from those shown in figure 9. 

Except where otherwise stated, the nominal translunar trajectory illustrated in 
figure 10 was employed throughout the analysis. The trajectory required 70 hours and 
37 minutes to reach perilune; it is plotted in a rotating-axis system to show the relative 
positions of the earth, moon, and vehicle at any given time. 

Preflight analy sis. - Again, it should be stated that the perturbations of r p and 
Vp from their nominal values are the principal contributors to the approach- guidance 
AV magnitude. Inasmuch as rp and Vp are functions of the deviation D (see fig. 3), 
these quantities can be used to determine the approach-guidance correction AV as a 
function of D. Determination of this relationship can be made before flight which means 
that the navigator need measure only the angles necessary to calculate D in order to 
determine the velocity correction required. 

The preflight procedure employed for determining AV as a function of D can be 
stated as follows: 

(1) By use of an n-body trajectory program, a number of perturbed trajectories (for 
example, 50) are generated from the point of the first midcourse correction to perilune, 
the perturbations being chosen randomly, based on the covariance matrix of errors after 
the first midcourse correction. During this process the second midcourse correction is 
applied at the aim point, no errors in execution being assumed. Actually, at midcourse 
the trajectories can be randomly perturbed about the nominal rather than around any or 
all trajectories perturbed by injection error; the effects would be the same in either case, 
as indicated in figure 3. Perturbing trajectories about the nominal eliminates the tedious 
task of generating a random sample of perturbed injection trajectories. This procedure 
also eliminates one step in the error analysis, that of having to account for the second 
midcourse maneuver. The small errors at second midcourse due to the perturbations at 
first midcourse would remain; however, their effect is negligible. 
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(2) For each such trajectory, the deviation D at the time of the approach-guidance 
maneuver is computed and rp and Vp are recorded. By using these data, rp and 
Vp are plotted as functions of D as in figure 3, and a curve is faired through the points. 
For onboard determination, D is calculated by equation (1) or (2); however, for preflight 
analysis, D can be calculated from the equivalent equation 

D = Z(x a - x n ) + m(y a - y n ) + n(z a - z n ) 

(3) For each value of D, corresponding values of rp and Vp are read from the 
faired curves. From these values of rp and Vp along with the corresponding values 
of r and V at the time of the approach guidance, a value of AV is computed by using 
equation (3) to yield a point AV,D. The resulting plot of AV as a function of D is 
used by the navigator to determine the required approach-guidance correction. 

Results of the foregoing preflight procedure are shown in figures 11 to 14 for sev- 
eral trajectories and times of approach-guidance measurement. Figure 11 presents the 
variation of perilune radius with the deviation D. For all cases hereinafter, the devi- 
ation is in the nominal orbital plane and in a direction perpendicular to the nominal 
range vector at the lunar sphere of influence. The deviations are determined for the 
times Tp in the figures. The first three parts of figure 11 pertain to the 70-hour trans- 
lunar trajectory illustrated in figure 10, the two different nominal values of rp being 
obtained by slightly changing the earth -injection conditions. Figure 11(d) applies to a 
90 -hour translunar trajectory. The slope of the curve in this figure is reversed because 
this trajectory is designed for counterclockwise motion about the moon. 

In the prediction of perilune radius by D, the scatter gives a good indication of the 
approach-guidance accuracy. Inspection of the plots for the various conditions shows that 
the scatter is not materially affected by the type of trajectory or the time of guidance 
measurement Tp. It can be noted in figure 11(c), however, that the scatter increases for 
perturbed trajectories falling below the lunar surface. Also, the variation in this figure 
is more nonlinear than in the other parts of figure 11; however, this nonlinearity has no 
effect on the guidance accuracy. The scatter in figure 11(d) is somewhat conservative in 
that it indicates more error than would ordinarily be obtained in practice. The extra 
scatter is due to the manner in which the midcourse-velocity errors were applied as dis- 
cussed in the section "Approach- Guidance Accuracy Characteristics." 

Figure 12 presents corresponding data on the variation of perilune velocity with the 
deviation D. It should be stated that the scatter for this quantity is not as important to 
the guidance accuracy as that for perilune radius. (Here again, in fig. 12(d), the slope is 
reversed because of the change in direction of the trajectory about the moon.) 
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Since flight-path angle appears in the preflight calculation of AV (eq. (3)), compar- 
ison was made between the two-body value of y as calculated from rp and Vp 
(eq. (4)) and the actual value determined directly in the n-body trajectory program. The 
results are shown in figure 13, where y and D are given for the time near the lunar 
sphere of influence (Tp = 14.617 hr). The two-body data are relatively smooth because 
they were determined from the faired rp and Vp curves in figures 11(a) and 12(a). 

The difference of about 0.7° between the two curves is due to two -body approximation. 

This difference makes it imperative to use equation (4) to calculate y for determining 
the AV, inasmuch as AV was derived from two-body theory. 

Figure 14 presents corresponding data on the variation with the deviation D of 
approach-guidance velocity as determined from the rp and Vp values in figures 11 
and 12. As previously stated, this variation is calculated by a preflight analysis of a 
number of perturbed trajectories; thus, the onboard computation of AV is eliminated. 

The signs for AV merely indicate whether the guidance-velocity pointing angle X is 
90° or -90°. Except for the angular measurements needed to calculate D, the variation 
of AV with D is the only information required by the navigator. For the trajectory 
of figure 14(d), the plot shows an error in AV of 0.3 m/sec at D = 0. The effect of this 
error, which may be due to the two-body approximation, can be corrected by merely 
offsetting the curve by 0.3 m/sec. 

The data in figure 14 are relatively smooth because faired values of rp and V p 
were used. The data of figure 15, however, were derived by using actual or unfaired 
values of r p (fig. 11(a)) and V p (fig. 12(a)) to determine AV; the scatter shown in this 
figure represents the velocity-correction error that would occur in practice because of the 
approximation. The deviation applies to a measurement time at the lunar sphere of 
influence (Tp = 14.617 hr); the dashed line, determined from data of figure 9, represents 
the AV requirement 5 hours after the time of measurement. 

Theoretical Considerations 

Error characteristics .- Position- and velocity-error ellipsoids representing char- 
acteristics of trajectory errors within the lunar sphere of influence due to onboard 
midcourse-guidance error are shown in figures 16 to 19. The guidance error is primar- 
ily due to a measurement error. Figure 16 gives the length of the axes of the position- 
error ellipsoid. The major axis is much greater than the other axes; hence, its direction 
may be considered to be representative of the direction of the position errors. The ori- 
entation of the major axis is shown in figure 17. In the upper plot, it is noted that the 
major axis lies always within 5° of the orbital plane. This feature provides the capability 
for determining the approach- guidance requirement by using only one star measurement. 
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The characteristics of the velocity-error ellipsoid are shown in figures 18 and 19. 
For figure 19, the major axis of the velocity-error ellipsoid is generally within 1° of the 
orbital plane. It can be seen in figure 19 that the error in the direction perpendicular to 
the velocity vector would be generally as large as the error along the velocity vector. It 
is the error in the perpendicular direction that has an effect on the perilune radius. As 
shown in figure 1, this velocity error can be several meters per second. Hence, it must 
be concluded that this error is correlated with D in order for the approach-guidance 
procedure to perform so well. 

Derivation of guidance relationships .- Attempts were made in the study to derive 
the relationship ip/D; however, no theoretical derivation was found. The only proce- 
dure shown thus far for determining the relationship between rp and D is by gener- 
ating a sample of trajectories randomly perturbed at first midcourse. This method 
requires calculations for a large number (for example, 50) of trajectories. In order to 
circumvent this problem, the following semiempirical procedure was developed for deter- 
mining the variation of rp with D; typical results are shown in figure 20. 

First, several magnitudes of D were arbitrarily selected, as shown by the sym- 
bols in figure 20. These values were then converted to magnitudes in the direction of the 
major axis of the position-error ellipsoid by dividing D by sin /3. (See fig. 17.) The 
direction cosines of the major axis were then used to obtain the change in the x, y, and 
z selenocentric position coordinates for the various points. The change in the vehicle 
velocity (error) for the various points was determined from the ratio of cr u /cr s in fig- 
ure 1. For example, for Tp of 14.617 hours, this ratio is 0.0074 m ^^ - c . The change 
in the nominal velocity vector was obtained by putting this error in the direction of the 
major axis of the velocity-error ellipsoid. (See fig. 19.) For the nominal trajectory 
used (fig. 10), which corresponds to clockwise motion about the moon, the selenocentric 
velocity magnitude should be reduced for positive values of D and increased for nega- 
tive values of D. This procedure would be reversed if the star direction were opposite 
to that in figure 2. 

For the semiempirical method, one needs the covariance matrix of midcourse 
errors and a knowledge of the pertinent error ellipsoids. This knowledge is obtained 
from trajectory-error -propagation programs by propagating the covariance matrix of the 
midcourse -guidance errors along the nominal trajectory. The orientation and shape of 
the error ellipsoids can be determined by simple matrix manipulation. (For example, 
see ref. 9.) The semiempirical method is only approximate because of the use of the 
major axes of the error ellipsoids to represent position- and velocity-error character- 
istics and because of the nonlinearity effects in propagating covariance matrices. The 
method, however, compares closely with the actual results, as shown in figure 20. 
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Appendix C contains several equations for calculating the ratio A(AV)/Arp. Fig- 
ure 21 presents actual data for this variation. Comparisons of actual and theoretical 
values given in figures 21(a), 21(b), and 21(d) show the equations to be very accurate. 

The theoretical expressions may be useful for various approach -guidance procedures in 
which rp is determined by some quantity other than D. (For example, in ref. 4, peri- 
lune radius is determined from measurements of the orbital angular velocity and body 
subtended angle.) It is of interest to note that the variation is extremely linear for maneu- 
vers made near the moon, (fig. 21(b)) and for trajectories designed to pass near the moon 
(fig. 21(c)). 

Combination of m aneuvers. - The foregoing figures correspond to the case where 
the approach-guidance maneuver is made after the second midcourse maneuver. Analysis 
of trajectories for which the aim point (time of second midcourse maneuver) was selected 
to be at nominal perilune time has shown that the approach-guidance procedure cannot be 
used without a second midcourse maneuver. Hence, the midcourse aim point must be 
chosen to be at or prior to the time of the approach -guidance measurements. The second 
midcourse maneuver can be conveniently combined with the approach-guidance maneuver, 
inasmuch as AVg is approximately linear with Ar mc and its direction is roughly the 
same for any injection error, as shown in figures 22 and 23, respectively. Each symbol 
represents a different perturbed trajectory due to injection error. Figure 22 represents 
position- and velocity- injection errors up to 10 km and 10 m/sec, respectively, and per- 
tains to an aim point at T = 56 hours (Tp = 14.617 hr). For these values, the figure 
shows that AVs can be predicted by Ar mc , the range measurement at first midcourse. 
For position errors much higher than 10 km the scatter would become excessive. Fig- 
ure 23 shows the precise in-plane angles of the second midcourse-guidance-velocity 
vector for the different perturbed trajectories. The angles are given with respect to the 
nominal velocity vector of the spacecraft and are shown to be roughly the same, especially 
at the larger values of Ar mc where the magnitude of AVg is relatively large. The 
degree of scatter shown in the figure is not significant; therefore, an average value can be 
used. The dispersions in the out-of-plane direction are even less than those shown in 
figure 23. 

The linear results of figure 24 were determined from each of the perturbed trajec- 
tories by applying the faired AVg values of figure 22 at an average angle of 15° from 
the nominal velocity vector. Figure 24, in effect, shows the change in perilune radius 
due to the second midcourse velocity AVg. Even though AVs would ordinarily be 
applied at angles ranging from about 10° to 35° (fig. 23), the small amount of scatter 
shown in figure 24 indicates that the perilune radius can be effectively corrected when 
AV is applied at a constant angle. Hence, the direction of AVs can be converted to 
a direction perpendicular to V n with little loss in accuracy. 
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To include the effect of AVg in the approach AV magnitude, which is in a per- 
pendicular direction, the following equation is used: 

AVadd = / Ar p \ /a(AV)\ 

Ar mc \^mc)\ Ar p / 

where Ar^Ar mc is obtained from figure 24 and A(AV)/Arp from figure 21(a). In this 
case 

^Yadd = -0.065 x 0.0200 = -0.0013 ™/sec 
ar mc km 

The quantity AVadd is an additional increment in the approach- guidance velocity and is 
added to V in the manner indicated in appendix B. In order for the navigator to include 
the effect of AVg in the approach-guidance maneuver, only a knowledge of the range 
measurement at first midcourse is required. It should be emphasized that if AVg is 
not applied in this manner, it must be applied prior to the approach guidance, inasmuch as 
no correlation exists between the deviation D and AVg. 

APPROACH-GUIDANCE ACCURACY CHARACTERISTICS 

In this section, the errors associated with the approach-guidance procedure are 
defined and analyzed. The analysis covers the range from the lunar sphere of influence 
to a point near the moon. It is assumed that there are no errors in the onboard calcula- 
tions due to human limitations. 

Effect of Approximation Error 

Table I summarizes the amount of perilune -magnitude error due to scatter for the 
various conditions illustrated in figures 11(a) to 11(d). This scatter error is an approx- 
imation error due to assumptions in the empirical approach-guidance procedure and is 
essentially the standard deviation of the difference between the data points and the faired 
line in each plot. The standard deviations given in the table, however, were obtained from 
a Monte Carlo analysis; that is, for each perturbed trajectory the faired value of AV 
was added to the corresponding V and the trajectory propagated to perilune. 

Conditions 1 and 2 differ only in the value of V used in calculating the guidance 
velocity (eq. (3)); as indicated, this difference has little effect on the guidance accuracy. 
Comparison of conditions 1 and 3 shows that the second midcourse correction can be 
combined with the approach-guidance maneuver with no loss in accuracy. Comparison of 
conditions 2 and 4 shows that the error due to scatter decreases as the approach-guidance 
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measurement is delayed to times closer to the moon. As the nominal perilune is lowered 
(condition 5), the scatter error increases. 

As previously stated, the results for condition 6 are conservative because of the 
manner in which the midcourse errors were simulated for this translunar trajectory. 

The errors were more nearly spherically distributed than would occur in practice from 
onboard midcourse measurements and lead to a higher inaccuracy in the approach- 
guidance procedure. In this case, the one-sigma values along the three axes of the error 
ellipsoid were 1.45, 1.11, and 0.168 m/sec. These values signify that the errors in the 
direction of the major and mean axes are roughly the same, and, in addition, the angle 
between the major axis and the spacecraft-velocity vector is about 83°. Figure 25 is an 
example of the approach-guidance scatter error for the 70-hour translunar trajectory 
with an approximate spherical distribution of velocity errors in the midcourse guidance 
application, which might occur in certain types of midcourse procedures. The error dis- 
tribution was obtained by incorporating errors at midcourse in each of the three compo- 
nents of the spacecraft nominal velocity. The errors were applied by including permu- 
tations of 0, ±1, and ±2 m/sec. For figure 25, CT r? p = 25.6 km, which, for all practical 
purposes, represents the absolute maximum scatter error for controlling approach tra- 
jectories by the present method. 

As a matter of interest, the scatter errors for some of the parameters at perilune 
are shown in table II. The off-nominal position and velocity deviations are cr s and op, 
respectively. The value oy p represents a change in the perilune-velocity magnitude; 
o u corresponds to changes in the three components of velocity (that is, to a change in the 
direction of the velocity vector). The reason for the large difference between these two 
values is evident in figure 19 which shows that the velocity errors at perilune are gener- 
ally in a direction 90° from the velocity vector. 

It is interesting to note that if perfect midcourse guidance is applied to the highly 
perturbed trajectory referred to in figure 3, the perilune-magnitude error is approxi- 
mately 25 km. This error is caused by the linear approximation made in using transition 
matrix theory. As indicated by one-sigma perilune errors listed in table I, the approach- 
guidance method corrects a large part of this error. 

Effect of Measurement Error 

One important effect on the approach-guidance accuracy is caused by the error in 
the required onboard measurements. The measurement -error equations are developed 
in appendix D. The nominal values of the measurement angles are given in figure 26 as 
a function of time to perilune. Also shown is the nominal variation of range to the moon. 
The change in the angle from the star to the moon results from the fact that the optimum 
star-measurement direction does not change as Tp decreases. 
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Figure 27 presents results on the effect of measurement error and the combined 
effect with approximation (scatter) error. The analysis pertains to a nominal 70-hour 
translunar trajectory with a perilune radius of 3403.6 km; however, the results would be 
about the same for other trip times and perilune radii. In the upper plot of figure 27, the 
solid curve shows the variation of a r , the one-sigma error in determining range by mea- 
suring the subtended angle of the moon. (See eq. (D3).) The one-sigma error in mea- 
suring the semisubtended angle was assumed to be 10 seconds of arc. The dashed curve 
shows the range-determination accuracy when the nominal value for range at the corre- 
sponding time is used; that is, the range is not measured. The error caused by using the 
nominal value is given by 

a r = a s cos j3 

where the position error a s is obtained from figure 1 and the angle /3 is given in fig- 
ure 17. As noted in the upper plot of figure 27, the dashed curve applies to trajectories 
with a fairly small midcourse-measurement error cr rjmc = 10 km. For u r rnc = 20 km, 
the values for the dashed curve would approximately double, and so on. The effect of the 
larger midcourse-measurement errors is shown by the other dashed curves in figure 27. 

The quantity is the one-sigma error in the deviation D due to the measure- 

ment errors cr r and og. The curves for oj) were obtained from equation (D2) by 
using values of cr r from the upper plot by assuming erg = 10 seconds of arc. As indi- 
cated in the lower plot of figure 27, the different dashed lines correspond to cases where 
r a - r n for three families of trajectories caused by different magnitudes of midcourse- 
measurement error. 

The data points shown for Arpy'AD were obtained from figure 5(b). This ratio is 
used to determine the effect of measurement error on the perilune-radius accuracy: 

Ar p 

a r,P = AD CT D 

The data points for the effect of scatter error were obtained from table I. 

The lower plot in figure 27 shows the perilune-radius error due to the combined 
effects of measurement and scatter error as determined by 
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It can be seen from the various plots that the effect of approach-guidance measurement 
error is small in comparison with that of scatter error. In fact, doubling the measure- 
ment error would not increase a rj p significantly. For the solid line in the lower plot of 
figure 27, the approach-range measurement is included, and, consequently, the magnitude 
of the midcourse -range -measurement error cr r>mc is insignificant. If the guidance mea- 
surements are made at distances relatively far from the moon, the range measurement 
does not substantially improve the approach-guidance accuracy, especially for low values 
of c r>mc . Actually, at the sphere of influence where the nominal value of 9 is 90°, the 
range -determination error has no effect on the approach-guidance accuracy. (See appen- 
dix D.) Also, if approach-guidance measurements are made closer to the moon, the nom- 
inal value of 9 at the sphere of influence may vary by several degrees with no increase 
in scatter error. (See fig. 6.) For example, from figure 26, it can be seen, that if 0 n 
at Tp = 14.617 hours were 88° or 92°, the time at which no range measurement would be 
required (0 n = 90°) would be shifted by about 4 hours. 

Effect of Maneuvering Error 

The effect of approach-guidance maneuvering errors on the perilune- radius accu- 
racy was examined. (See figs. 28 and 29.) Since maneuver timing error represents 
error in the direction of range, it is considered to be negligible, especially at or near the 
sphere of influence where the range error has no effect. In fact, timing error need not be 
considered in any type of onboard guidance procedure inasmuch as errors of 1 or 2 min- 
utes applied to the midcourse procedure of reference 9 lead to aim -point errors of only 
several kilometers. 

Pointing error .- The effect of pointing-direction error is shown in figure 28. No 
statistical analysis was performed; rather, the individual effects of the in-plane and out- 
of -plane components of this error were determined. The plane referred to is the instan- 
taneous nominal earth-moon-vehicle plane which is essentially the selenocentric orbital 
plane of the vehicle. The data were calculated for only one perturbed trajectory, but are 
representative of any perturbed trajectory. The results indicate that the effect of pointing 
error can be considered to be negligible; that is, errors of several degrees would have no 
great effect on the approach-guidance accuracy when the nominal X is selected as 90°. 

Velocity-cutoff error .- Inasmuch as the approach-guidance AV magnitude is rel- 
atively low, especially near the lunar sphere of influence (fig. 15), the effect of guidance- 
velocity-cutoff error was investigated statistically and is shown in figure 29. The 
guidance- velocity- requirement curve in the upper plot used in determining this effect is 
taken from figure 9. 

The curve showing the effect of cutoff error on the perilune -radius error was deter- 
mined from the inverse of the velocity requirement by the following equation: 
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a 


r,P 


Ar g - (o,2 x 10' 3 ) 

A(AV) V 


(5) 


where the value 0.2 is the one-sigma error, in meters per second, assumed in cutoff 
velocity. 

The lower plot in figure 29 shows the variation in the overall error in perilune 
radius with time of terminal-guidance measurement, as determined by the equation 


CT r,p 


2 

( CT r,p) cutoff + 



2 

scatter 


+ 


/ Ar p 

UD°b, 



In general, 40 percent of the error is attributed to each of the first two sources and the 
remaining 20 percent to the measurement error. The cutoff error in figure 29 applies 
when X = 90°. Application of the AV vector at other values of X (fig. 8) reduces the 
effect of cutoff error as is shown by equation (5). Another error in the guidance -velocity 
magnitude, called the proportional error, is a small constant percentage of AV and was 
found to be negligible. 

The effect of a cutoff error of 0.2 m/sec was also determined for some of the other 
quantities at perilune for condition 2 given in table n. The effect was to increase oy p 
from L99 to 2.90 m/sec and to increase a s and o u to only 27.20 km and 11.40 m/sec, 
respectively. 

With respect to times near the sphere of influence, the perilune-radius error can be 
reduced by about one-half if approach guidance is delayed to 5 hours from the moon. The 
fuel requirements, however, will be tripled, as shown by the upper plot in figure 29. 


Effect of Ephemeris Error 

Even though the lunar ephemeris error is very small (1 or 2 km), its effect on the 
approach-guidance accuracy was examined. Larger ephemeris errors were included to 
give some indication of the effect on the approach-guidance procedure when controlling 
interplanetary trajectories. Ephemeris error affects the accuracy of the guidance pro- 
cedure because the guidance measurements are referenced to a nominal trajectory which, 
in turn, is based on a certain location of the moon. This type of error was introduced 
into the procedure by changing the ephemeris of the moon in the n-body trajectory pro- 
gram. In figure 30, it is seen that the approach-guidance procedure compensates, to a 
large extent, for the effect of ephemeris error. For example, for ephemeris error in the 
direction of the moon's motion, the error in controlling perilune radius to the nominal 
value is only about 20 percent of the ephemeris error. The angle between the nominal 
range vector and the earth-moon line is about 51°; therefore, the perilune error would be 
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different if the ephemeris error were in the direction perpendicular to the nominal range 
vector. In fact, this direction generally has the greatest effect on interplanetary 
approach -guidance error, even though it is the direction for minimum ephemeris error 
(ref. 10). Although trajectory characteristics near a planet are not identical to those 
within the lunar sphere of influence, the results in figure 30 would roughly apply to inter- 
planetary flight. 


Prediction of Perilune Time 

The variation of time of perilune passage with approach-guidance velocity is shown 
in figure 31 for a number of perturbed trajectories. As stated in the figure, the perturbed 
trajectories are due to one-sigma midcourse- range-measurement errors of 22 km. The 
magnitude of this error, however, affects the magnitude of approach AV and not the 
degree of scatter in the data. A prediction of this time, based on the magnitude of 
approach AV, may be important for deboost into lunar orbit. The scatter of the data 
about the faired line indicates a one-sigma prediction capability of about 82 seconds. The 
maximum error shown is only about 180 seconds. 

CONCLUSIONS 

A method for determination of approach-velocity corrections for a space vehicle 
has been presented. The method has been applied to the approach phase of earth-moon 
trajectories; however, it will also apply to reentry control for moon-earth trajectories. 

The method is unique in that only a single onboard position measurement is required to 
determine the guidance correction for controlling the perilune -radius magnitude. This 
position measurement is in the form of a deviation from a nominal trajectory and normally 
requires a subtended-angle measurement and a star-to-body angular measurement. Use 
is made of preflight calculations of the nominal trajectory and of various parameters for 
trajectories perturbed about this trajectory. These calculations not only provide the 
AV-magnitude variation with deviation, but preselect the measurement star and direction 
of the velocity-correction vector. These preflight determinations are the only require- 
ments for the approach method, other than performing the measurements. The charac- 
teristics of the method are such that perilune accuracy of around 10 km can be obtained 
even though large differences may exist between the actual and nominal trajectories. 

From the error analysis performed on the method, the important results concerning the 
onboard approach-guidance method are: 

1. The method is especially applicable to controlling errors resulting from onboard 
midcourse -measurement procedures. 

2. The method can be applied anywhere within the lunar sphere of influence. 
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3. The measurement star must be in a specified direction for adequate guidance 
accuracy. This direction is perpendicular (±2°) to the direction of the nominal range 
vector at the lunar sphere of influence and does not change with time within the sphere. 
For highest sensitivity, the star should be near the orbital plane; however, angles up to 
30° or more could be tolerated. 

4. A range determination (for example, by measuring the subtended angle of moon) 
is not necessary if the guidance maneuver is made far from the moon, that is, near the 
lunar sphere of influence. 

5. For one-sigma errors of 10 seconds of arc in the onboard angular measurements, 
the one-sigma error in controlling perilune radius is approximately 13 km. This error 
includes a one-sigma error of 0.2 m/sec in the cutoff velocity and corresponds to times 
near the lunar sphere of influence. If approach guidance is delayed to 5 hours from the 
moon, the perilune error is reduced by about one -half. The fuel requirements, however, 
are tripled. 

6. Only a small part of the perilune-radius error is attributed to measurement 
error, if it is assumed that the one-sigma measurement error is 10 seconds of arc. In 
general, the cutoff-velocity error contributes 40 percent of the error, as does the effect 
of approximation error associated with the empirical procedure. 

1 . The approach-guidance procedure is insensitive to other types of maneuvering 
errors such as timing and pointing direction. 

8. Whereas the method is designed principally to control the magnitude of the peri- 
lune radius, it has been shown also to control perilune position and velocity within reason- 
able limits. 

9. It has been shown that the present guidance method compensates, to a large 
degree, for the effect of ephemeris error. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 24, 1970. 
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APPENDIX A 


EQUATIONS FOR UPDATING GUIDANCE MEASUREMENTS 

The two angular measurements, 9 and a, normally used in determining the 
approach-guidance velocity must be referenced to a common time. With the type of mea- 
surement scheme shown in the sketch, both measurements could be made essentially at 



the same time so that updating would not be necessary. However, if some other measure- 
ment scheme were employed or if more than one set of measurements were needed to 
reduce the measurement error, updating the measurements could be required. 

Either measurement can be updated to a common time (small increments of time 
only) by means of data such as that presented in figure 32 and by the following equations: 


O-SmH-f At 


a = Qf m + — At 
dt 


where At is the increment between the measurement time and the common time. In con- 
verting a measurement to a common time, the corresponding value of the rate of change of 
angle with time along the nominal trajectory, d0/dt or da/dt, must be selected for the 
given time T p . The value of T p would need to be known only approximately; hence, the 
nominal value would suffice. The values of d0/dt and da/dt, shown in figure 32, are 
precalculated from the equations which follow. 
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APPENDIX A — Continued 


The rate of change of the measurement 0 (angle between star and moon center) 
with time along a nominal trajectory is given by 

d0 90 dx 90 dy [ 90 dz 
dt 9x dt 3y dt 9z dt 

where dx/dt, dy/dt, and dz/dt are known values for the nominal trajectory in the 
selenocentric coordinate system and, as can be determined from equations given in 
reference 11, the partials for the angular measurements 0 are 


9 8 _ 

x(my + nz) -£(y 2 + z 2 ) 

dx 

r 2 j~r 2 - (lx + my + nz)2^j l/ 2 

3 9 

y (lx + nz) - m(x 2 + z 2 ) 

9y 

r 2j^ r 2 - (lx + my + nz)2] 

90 _ 

z (lx + my) - n(x 2 + y2) 


9z r 2 £r2 - (lx + my + nz PJI72 


If the star is in the orbital plane, the equation for d0/dt is equivalent to the rate of 
change of true anomaly with respect to time. 

The rate of change of the measurement a (semisubtended angle of moon) with time 
along a nominal trajectory is given by 

da?_ 9a; dx ( 9a dy t da dz 
dt 9x dt 9y dt 9z dt 

where, again, the derivatives with respect to time are commonly known values for the 
nominal trajectory and the partials for the angular measurement a, which can be deter- 
mined from the relation 

a = tan“l ^ 

/r 2 - R 2 


are 


da 

dx 


-X 


R 

r 2 \/r2 - r2 


25 



APPENDIX A - Concluded 



26 



APPENDIX B 


DERIVATION OF VELOCITY CORRECTION REQUIRED TO 
ATTAIN DESIRED PERILUNE 

From the law of conservation of angular momentum 

rp,n 2 V P)I1 2 = r2v^ 2 cos 2 y^ = t 2 Va 2 cos 2 (y + Ay) 

where y is the flight-path angle prior to the approach-guidance correction, and y A 
and V A represent values immediately following the correction. 

Expanding equation (Bl) gives 

r p,n 2 Vp,n 2 = r 2 VA 2 (cos y cos Ay - sin y sin Ay) 2 



Solving for Ay from the sketch gives 


2i T 9 2i r 9 cos y(V + AV cos A) sin y AV sin A 

V n n 4 * ~ I* V * __ ~ t t 


Pj n *P,n 


rp jn 2Vp >n 2 = r 2 [cos y(V + AV cos A) - sin y AV sin a] 


Expansion yields 


rp )n 2Vp jn 2 = r 2 [cos 2 y(V + AV cos A) 2 + AV 2 sin 2 y sin 2 A 
-2 AV(V + AV cos A) sin y cos y sin xj 


r p n 2 Vp jn 2 = r 2 [AV 2 cos 2 (y + A) + 2V AV cos y cos(y + A) + V 2 cos 2 y] 
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APPENDIX B - Continued 


But 


V 2 _ 2 M - y 2 _ 
V P’ n r p,n A r 


(B4) 


and from the sketch, it is seen that 

Va 2 = V 2 + AV 2 + 2V AV cos X 


(B5) 


Substituting equations (B4) and (B5) into equation (B3) yields 


r 2 (v 2 + AV 2 + 2VAV cos\-^ + ^ 
P) n \ r rp >n/ 


= AV 2 r 2 cos 2 (y + X) + 2V AVr 2 cos y cos (y + X) + V 2 r 2 cos 2 y 


whence 


AV 2 


r D n 2 - r 2 cos 2 (y + X) 


p,n 

+ V 2 fr p5 n 2 - r 2 cos 2 y ) + r Pj n 2 (^ " ¥) = 0 


+ 2V AV 


r p n 2 cos X - r 2 cos y cos(y + X) 


Solving for the velocity correction yields 


AV = 


V £r 2 cos y cos(y + X) - r Pjn 2 cos xj 


rp,n 2 “ r 2 cos 2 (y + X) 


r p,n{ ( r2 " r p,n 2 ^ 2 sin 2 X + 


r 2 cos 2 (y + X) - r p>n 2 


2M 2 M V|1/ 2 


>,n" rJJ 


rp n 2 - r 2 cos 2 (y + X) 


(B6) 


The value of y is derived from 


r P V P 
cos y = -^r- 


(B7) 


where r p and V p are the predicted perilune values. (See figs. 11 and 12.) The alter- 
nate signs of the second term in equation (B6) correspond to correcting to either side of 
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APPENDIX B - Concluded 


the moon. The sign which results in the lesser value of AV would ordinarily be chosen. 
It should be noted that AV is added to (subtracted from) V if the nominal trajectory is 
designed to rotate clockwise (counterclockwise) about the moon (viewed from northerly 
direction) as shown in figure 2. 
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APPENDIX C 


VELOCITY REQUIREMENT FOR APPROACH GUIDANCE 

Several analytical expressions for the variation of approach-guidance velocity AV 
with perilune radius rp are presented in this appendix. The velocity requirement 
A(AV)/Arp is useful for error analysis (see fig. 29) and as mentioned in the main text, 
for obtaining AV if r p can be determined directly from some type of approach- 
guidance measurement. Values resulting from each of the following expressions were 
essentially the same and are compared with actual values in figure 21. 

The first equation presented is the exact expression for velocity requirement (at 
small values of AV) and was determined by differentiating equation (B6) with respect to 
rp. The resulting equation is 


2r2Vrp cos(y + X)Q:os y + cos X cos(y +xfj 
® rD fr2 cos2(y + X) - rp2j 2 


(rV sin 


in X)2 - (rpV sin X)2 + Jr2 cos2(y + X) - r p^j^y- - "r“Jj cos ^(> / + X) + rp 2 J 

jr 2 cos 2 (y + X) - Tp^j 

r p3v2r sin2x + |U r 2 cos2(y + X) + jump 2 - 2/urp^ 
p[r2 cos2(y + X) - r p 2j |(rV sin X)2 - (r p V sin X) 2 + [r 2 eos 2 (y + X) - Tp^J 


Simplified expressions which were derived from work accomplished in reference 12 


are 


9V 

9r„ 


= ± 


V p 6 

1 +e 


1 - 3 + e 

fa) 

2 o / 
1 + - 

' r p\ 3 

L 1 +e' 

\ r / 

1 + e\ 



1/2 


for optimum-angle thrust direction and 

V P e 


av 


3r 


P 


± 


1+6 



r 


for thrust in the direction of the local horizontal (which is near optimum). 


30 



APPENDIX D 


MEASUREMENT-ERROR EQUATIONS 

Derived in this appendix are the equations used in analyzing the effect of measure- 
ment error in the approach-guidance procedure. The measurement equation is 

D = r m cos 0 m - r n cos 0n = r m cos 0 m - Constant 

where D is the position deviation calculated from the measurements r m and 0 m . It 
is seen that 

dD = cos 8 dr - r sin 9 d0 (Dl) 


or 


AD = cos 9 Ar - r sin 8 A 8 

Uncorrelated Measurements 

If it is assumed that Ar and A 9 are random uncorrelated errors, 

ct d 2 = (cos 9 cr r ) 2 + (r sin 9 cfq)2 (D2) 

where nominal values for r and 9 can be used. The quantity oq is normally taken 
as constant, whereas cr r is a function of range. For times near Tp = 14.614 hours 
(0 ® 90°), 

crj-)2 ~ r 2<jg2 

where cr g is the standard deviation of the star-to-moon angular measurement error. 
Inasmuch as the range -measurement error is insignificant for 9 ~ 90°, the nominal value 
can be used in place of the range measurement when D is calculated. 

For times closer to the moon, the range measurement, and hence its error a r , 
become significant. The error in range determination is caused by error in the 
subtended -angle measurement a and uncertainty in the knowledge of the moon's 
radius R. From the equation 

R = r sin a 
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APPENDIX D — Continued 


it is seen that 

dR = r cos a dor + sin a dr 


or 


sin a Ar = AR - r cos a Aa 


so that 



Now, since a is small 


or 


^ = (s“R) 2 + ( r H 2 

«r 2 = ^ 2 ~V) 


(D3) 


Correlated Measurements 

If the angles 9 and a are measured by the method suggested in appendix A, their 
errors would be correlated. From the equations given in the sketch in appendix A, it is 
seen that 


d 9 = |(d0i + d02) 
da; = “ d02j 


Also 


r 


R 

o? 


so that 


dr= -^ da = ||(d« 2 -dei) 
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APPENDIX D - Concluded 


If the relations for dr and d0 are substituted into equation (Dl), then 

co - (xK" 2 ' d 9 i) - !L f Ll ( d91 + d82 ) 


dD = 


or 


AD = 


' r cos 0 - sin 6] A0 2 - (*-££±1 + sin 0) A0j. 


If it is assumed that 0j and 0 2 are measured separately, that is, A0j and A02 are 
random uncorrelated errors with equal variances 


o 2 _ ££ 

a D - 4 


2 cos 2 0 + sin 2 #) og 2 


(D4) 


or 

a D = ^ r^sin 2 0 + ^ cos 2 ©) Oq 


For times near Tp = 14.614 hours (0 ~ 90°), 


Opj ~ 0. 707rog 


It is interesting to note that in this case, the errors are smaller than those for the previ- 
ous case (eq. (D2)). 

In equation (D4), the error in R has not been taken into account. Including the 
error in R gives 


ct D “ 



r 2 2 \ 1/2 
+ R 2 Gr ) 


where o D is determined from equation (D4) and the prime denotes the inclusion of the 
uncertainty in lunar radius. 
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TABLE I. - EFFECT OF SCATTER ERROR ON ACCURACY OF 
CONTROLLING PERILUNE DISTANCE 


Condition 

Nominal 
trip time, 
hr 

^P, n » 

km 

T 

A p,n> 

hr 

Value of V 
used in equation (3) 

Second midcourse- 
approach- guidanc e 
maneuvers 

°p,n> 

km 

1 

70.617 

3403.6 

14.617 

v n 

Separate 

8.78 

2 

70.617 

3403.6 

14.617 

v a 

Separate 

8.60 

3 

70.617 

3403.6 

14.617 

V n 

Combined 

8.27 

4 

70.617 

3403.6 

4.617 

v a 

Separate 

5.5 

5 

71.07 

1892.8 

15.07 

v n 

Separate 

10.52 

6 

90 

2237.3 

15.0 

v n 

Separate 

14.79 


TABLE II.- EFFECT OF SCATTER ERROR ON ACCURACY OF 
CONTROLLING PARAMETERS AT PERILUNE 

[Condition 2 of table i] 


Parameter ° 

Inclination, deg 0.36 

Longitude, deg 0.37 

Latitude, deg 0.17 

Vp , m/sec 1.99 

s, km 25.32 

u, m/sec 11.14 
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Perilune radius, r„, km Perilune radius, 



(a) Star in direction of 
selenocentric velocity 
vector (close to earth- 
moon- vehicle plane). 


(b) Star in direction of 
major axis of position- 
error ellipsoid (about 
^-7° out of earth- moon- 
vehicle plane. 



D, km 

(c) Star in earth- moon- vehicle 
plane and in direction per- 
pendicular to hyperbolic 
asymptote (about 3.5° and 9-5° 
from perpendicular to seleno- 
centric velocity and radius 
vector, respectively). 



(d) Star in earth- moon- vehicle 
plane and in direction per- 
pendicular to selenocentric 
radius vector. 


Figure 4.- Comparison of scatter at Tp = 1^.617 hr for deviations 
in different directions. cr r mc = 10 Km. (Lines are faired 
through data points.) 



Perilune radius, 







Percent increase in AV 








1 800 


Lunar radius 



D, km 

(c) 70 -hr translunar trajectory; r p,n = 1892-8 km; Tp = 15-07 hr 


Figure 11. - Continued. 





Perilune velocity, 



(a) 70 -hr translunar trajectory; r p,n = 3 ^ 03.6 km; 
V p ^ n = 1.9805 km/sec; T p = 14-.617 hr. 

Figure 12.- Variation of perilune velocity with D for 
50 perturbed, trajectories. a T mc = 22 km. 





Peri ui 



Figure 12.- Continued. 





Calculated flight-path angle, y, deg 













300 -600 -400 -200 0 200 


D, km 

Figure 15*- Variation of approach-guidance-velocity requirements 

with D for two times. Measurements made at sphere of influence 



Length of axis, km 



Time to nominal perilune time, T hr 

Figure l6.- Magnitude of axes of one-sigma position-error ellipsoid 
resulting from measurement error in onboard midcourse-guidance 
procedure. cr r mc = 10 km; ct 0 = 10 seconds of arc. 


In - plane angle, (3, deg 


d cn u 


CO CD 




Figure 17. - Orientation of major axis of position-error ellipsoid caused 
by measurement error in onboard midcourse- guidance procedure. 




Time to nominal perilune time, T hr 

Figure 18.- Dimensions showing shape of one-sigma 

velocity-error ellipsoid resulting from measurement 
error in onboard midcourse-guidance procedure. 

°r,mc = 10 km; tfg = 10 seconds of arc. 




Angle between major axis and velocity vector, deg 






-600 -400 -200 0 200 400 600 

D, km 


Figure 20.- Comparison of actual variation and calculated variation 
of r_ with D. T-. = 14.617 hr. 


Guidance velocity, A V, m/sec 



4400 


(a) 70 -hr translunar trajectory; 

r p,n = 3^03.6 km; T p = l>i.6l7 hr. 

Figure 21.- Variation of approach-guidance velocity 
with perilune radius for 50 perturbed trajec- 
tories. o T>mc = 22 km. 




Figure 21.- Continued. 









Midcourse velocity correction, AV> m/sec 





Figure 23-- Variation in direction of second midcourse-velocity vector. IL-p = 9-5 hr. 




Perilune radius, r n , km 



Figure 25-- Prediction of perilune radius resulting from approximately spherical 
distribution of midcourse- velocity errors. Maximum velocity error, ±2 m/sec 
in any component; D at 


T p = 12.117 br. 


Range to moon center, r, km 
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Time to nominal peri lune time, T p 

Figure 26.- Characteristics of nominal trajectory. Star is 
plane and in direction perpendicular to range vector at ' 
influence. 




A . Range- determination 

scatter’ km p/ AD °D' km error, o km 




Time to nominal perilune time, T hr 

Figure 27-- Approach-guidance-accuracy characteristics due to measurement error and 
scatter error. Measurement errors a a and 00 = 10 seconds of arc; moon-radius 
uncertainty, 0.8 km. Range to moon measurement used for solid lines in upper 
two plots and lowest plot. 
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Figure 29- - Approach-guidance error with effect of velocity-cutoff error included. 
One-sigma value of velocity-cutoff error, 0.2 m/sec. Range to moon measurement 
used for solid line in bottom plot. 


Perilune-radius error, Ar n , km 



Figure 30.- Effect of lunar ephemeris error on approach-guidance 
accuracy. Tp = 12.117 hr. (For ephemeris error in direction 
opposite to moon's motion, direction of Arp is away from moon.) 
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Time to nominal perilune time, 



hr 


Figure J2.- Rate of change of measurement angles with time 
along nom'nal trajectory. 
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